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Abstract 



Within the framework of the dinuclear system model, the production of superheavy ele- 
ment Z=117 in possible projectile-target combinations is analyzed systematically. The calcu- 
lated results show that the production cross sections are strongly dependent on the reaction 
^ ' systems. Optimal combinations, corresponding excitation energies and evaporation channels 

are proposed in this letter, such as the isotopes in ^^Ca induced reactions in 3n 

evaporation channels and the reactions ^^^c +246,248 ^^^^ channels, and the 

system V -|-244 • channel. 

PACS: 25.70.Jj, 24.10.-i, 25.60.Pj 

The synthesis of superheavy elements has obtained much progress experimentally using cold 
fusion reactions with double magic nucleus ^osp^ or nearly magic nucleus '^^'^Bi as targets and hot 
fusion reactions with double magic nucleus ^^Ca bombarding actinide nuclei in fusion-evaporation 
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reaction mechanism. ^ Especially in recent years, superheavy elements Z=113-116, 118 have 
been synthesized successfully in laboratories.^^ Of course, further confirmations on these new 
elements have to been done. Various theoretical models have also been proposed to describe 
the formation of superheavy nucleus.^ ^ ^ ^ ^ ^ Among these models, the dinuclear system 
(DNS) model can describe very well a series of available experimental data.^^ We have further 
developed the DNS model by taking the diffusion process to be coupled with the relative motion, 
and by introducing barrier distribution function method in the capture of two colliding nuclei and 
in the fusion process.^ The decay of the formed DNS (quasi-fission) and the fission of heavy 
fragment are also considered in the model. For all of reaction systems including cold fusion, 
^^Ca induced reactions and other projectile-target combinations etc, we can use the same fixed 
parameters to reproduce experimental data. 

In the DNS concept, the compound nucleus is formed by nucleon transfer from the light nucleus 
to the heavy one by overcoming the inner fusion barrier. The evaporation residue cross section 
can be written as a sum over all partial waves J at the centre-of-mass energy Ecm.i 

Jmax 

(^ER{Ec.m) = 7^—^ X! (^"^ + ^)T{Ec.m.-, J)PcN{Ec.m., J)Wsur{Ec.m. , J)- (1) 

Here T{Ec.m., J) is the probability of two colliding nuclei overcoming the potential barrier in the 
entrance channel to form the DNS. PcN{Ec.m.i J) is the probability that the system will evolve 
from a touching configuration into the formation of compound nucleus against quasi-fission of the 
DNS and fission of the heavy fragment. The last term is the survival probability of the formed 
compound nucleus, in which the fissile compound nucleus will decay by emitting 7 rays, neutrons 
and charged particles etc. In the calculation we take the maximal angular momentum as Jmax = 30 
since the fission barrier of the heavy nucleus disappears at very high spin.E^ 

The same as in the nucleon collectivization model, ^ the transmission probability is calculated 
by introducing a barrier distribution function, which can reproduce very well available experimen- 
tal capture cross sections.^ ^ The fusion probability is obtained by solving Master Equation 
numerically in the potential energy surface of the DNS. The time evolution of the distribution 
function P{Ai^e\,t) for fragment 1 with mass number Ai and local excitation energy e\ is de- 
scribed by the following Master Equation, 

'^^^^'f'^^ = Y: Wa„a', it) [dA,PiA[,er, t) - dA'P{A,,el, t) - 
at 

A''^(e(t)) + A^^^(0(t))]p(Ai,el,t). (2) 

Here Wai,a[ is the mean transition probability from the channel {Ai,el) to {A[,el'), while (Iai 



denotes the microscopic dimension corresponding to the macroscopic state {Ai,el). The sum is 
related to Wa-^^^a'^ and taken as = Ai ± 1 for one nucleon transfer. The local excitation energy 
e\ with respect to fragment Ai is related to the intrinsic excitation energy of the composite system 
and the driving potential of the DNS. The intrinsic excitation energy is provided by the kinetic 
energy of the relative motion, in which the barrier distribution of colliding system is introduced 
in the dissipation process. The motion of nucleons in the interacting potential is governed by 
the single-particle Hamiltonian.fl^ Quasi- fission rate A^-'^ and fission rate (for heavy fragment) 
are estimated with one dimensional Kramers formula. [^^^ After reaching the interaction 
time in the evolution of P{Ai, e\,t), the formation probability at Coulomb barrier B and angular 
momentum J is given by 

Abg 

PcN{Ec.m.,J,B)= P{A,,el,nnt{Ec.m.,J,B)), (3) 

Ai=l 



where the interaction time Tint is obtained by deflection function method. ESI So we can obtain 
fusion probability as 

Pcn{Ec.^., J)= J fiB)PcN{Ec.m., J, B)dB. (4) 

For the barrier distribution function f{B), we take an asymmetric Gaussian form in the calculation. 
The survival probability of the thermal compound nucleus is given by 

where the E^^, J are the excitation energy and spin of the compound nucleus, respectively. The 
realization probability P{E^^,x, J), neutron evaporation width and fission width Fj are cal- 
culated with the statistical evaporation theory. 

The partial cross sections of the capture, fusion and evaporation residue are written respectively 

as 

<Jcap{Ec.m., J) = 2^^(2J + l)T{E,,„,,, J), (6) 
C^fusiEc.m., J) = 2^E^ra ^ ^)T{Ec.m., J)PcNiEc.m., J), (7) 

aevr{Ec.m., J) = 2j^^{'^J + ^)T{Ec.m., J) PcN{Ec.m. , J)Wsur{Ec.m. , J)- (8) 

Within the framework of the DNS model, we have analyzed the partial cross sections in the 
three stages for the reaction '^^Ca +'^'^^ Cm at incident cm. energies 201.91 MeV (solid line) 
and 211.91 MeV (dashed line) as shown in Fig.l (a)-(c), corresponding to excitation energies of 
compound nucleus 35 MeV and 45 MeV respectively with the relation E^^ = Ec.m. + Q- The 



reaction Q value is given by Q = AMp + AM^ — AMq, and the corresponding mass defects 
are taken from Ref.[T5j for projectile, target and compound nucleus, respectively. We can see 
that the dependence of the partial cross sections on angular momentum is quite different. At 
the considered excitation energies, higher angular momentum has larger contribution for capture 
partial cross section. However, the maximal positions of the fusion and evaporation residue partial 
cross sections begin to remove towards lower angular momentum since fusion probability and 
survival probability decrease rapidly with increasing angular momentum. In Fig.l (d) we take 
a sum of the evaporation residue partial cross section over maximal angular momentum until 
Jmax = 30 in 2n-4n evaporation channels. Calculated results can reproduce experimental data 
except for the 4n evaporation channel which predicts too high excitation energy. It may be 
related to the incident energy loss in the target or to the calculated Q value. As a whole, the 
maximal values of the calculated aER are consistent with the experimental results. 

With the same method as in Fig.l (d), we have calculated the evaporation residue excitation 
functions of the cold fusion reaction ^"^Se +^°^ Bi and the reaction ^^Ca +^^^ Bk in 2n-4n channels 
to synthesize superheavy element 117 as shown in Fig. 2. The maximal production cross sections 
are 0.015 pb and 0.86 pb for the reactions ®^S'e +'^^^ Bi in In channel and ^^Ca +^^'' Bk in 3n 
channel, corresponding to excitation energies 12 MeV and 37 MeV, respectively. 2n channel and 
4n channel have smaller production cross sections, which are 0.61 pb and 0.48 pb, respectively, 
corresponding to excitation energies 26 MeV and 47 MeV. Therefore, it is favorable to synthesize 
new superheavy element 117 in 3n evaporation channel for the systems ^^Ca +^^^ Bk. Isotopic 
trends to produce superheavy element 117 in cold fusion reactions based target and in 

^^Ca induced reactions are analyzed using the DNS model in Fig.3. It is shown that isotope 
projectiles ^^5*6, "^^Se and isotope targets '^'^^Bk, '^'^^Bk have larger production cross sections. ^^Ca 
induced reactions have larger production cross sections than cold fusion reactions. Here we only 
pay attention to 3n evaporation channel since it has larger cross section than other evaporation 
channel. Isotopic dependence of production cross section are mainly affected by driving potential, 
neutron separation energy and shell correction energy. Usually the incident channel with more 
asymmetric combination (larger mass asymmetry degree of freedom) has larger fusion cross section 
due to smaller inner fusion barrier. Smaller neutron separation energy and higher fission barrier 
(mainly coming from shell correction for superheavy nucleus) will increase the survival probability 
of the formed compound nucleus. The reaction ^^Se +^°^ Bi forms the element 117 with 172 
neutrons, which is a sub-magic number. In the case the one neutron separation energy is larger, 
which leads to smaller survival probability, so as to smaller evaporation residue cross section. 
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Neutron-rich combinations will reduce neutron separation energy, but it does not guarantee the 
fusion probability to be larger. There exists the competition between Pcn and Wsur as a function 
of the neutron richness of the system. Pairing effect also affects the neutron evaporation width 
and the fission width. It may be good way to synthesize superheavy nucleus using radioactive 
nuclear beam in the future due to larger survival probability of formed compound nucleus. 

Superheavy elements have been successfully synthesized using fusion-evaporation reaction in 
laboratories in the world. Usually double magic or nearly double magic nuclei are selected to 
synthesize superheavy nucleus because these combined systems have larger Q value, which can 
reduce excitation energy of compound nucleus so that fusion probability and survival probability 
increase, such as ^^Ca induced reactions in Dubna, and cold fusion reactions based on ^°^P6 or 
'^^^Bi target at GSI. We calculated other reaction systems to produce superheavy element 117 
as shown in Fig. 4, such as (a) ^^Sc+^^^ Cm(sohd hne) and '^^Sc-\-^'^^ Cm(dashed line), as well 
as (b) +2^^ Pu, (c) 55 Mn +'^^^ U and (d) ^'^Co +^32 j,^ -^^ 2n-4n channels. These nuclei 
have longer half-lives and are possible to use in laboratory with intense beams. Calculations give 
that the reactions '^^Sc -l-^^^ Cm and '^^Sc -l-^^^ Cm have larger production cross sections in 3n 
and 4n evaporation channels than other systems since lower inner fusion barriers can increase 
the fusion probability, which are suitable to synthesize the new element. However, the target 
material ^^^'^"^^Cm are very difficult to be prepared. To prepare '^'^'^Pu is relatively easy, but one 
has to face the smaller production cross section of 0.6 pb. In Table 1 we list optimal excitation 
energies, maximal evaporation residue cross sections as well as transmission, fusion and survival 
probabilities at J = in ln-4n evaporation channels for all of the possible combinations to produce 
the element 117. 

In summary, possible combinations are analyzed to synthesize new superheavy element 117 
within the framework of the DNS model. The partial cross sections in the capture, fusion and 
evaporation process for the reaction ^^Ca -l-^^^ Cm are calculated as a test of the model. Within 
error bars experimental data can be reproduced very well. Isotopic trends of the production 
of clement 117 in cold fusion reactions based ■^^^'^Bi target and in "^^Ca induced reactions are 
investigated systematically. Calculated results show that the reactions ^^Ca +248,249 
evaporation channels and '^^Sc +246,248 channels are favorable to produce the 

element. The corresponding excitation energies are also given in the letter. 



References 

[1] Hofmann S, Miinzenberg G 2000 Rev. Mod. Phys. 72 733 

[2] Oganessian Yu Ts, Yeremin A V, Popeko A G et al 1999 Nature 400 242 

Oganessian Yu Ts, Utyonkov V K, Lobanov Yu V et al 2000 Phys. Rev. C 62 R041604 
Oganessian Yu Ts, Utyonkov V K, Lobanov Yu V et al 2001 Phys. Rev. C 63 R011301 
Oganessian Yu Ts, Utyonkov V K, Lobanov Yu V et al 2004 Phys. Rev. C 69 R021601 
Oganessian Yu Ts, Utyonkov V K, Lobanov Yu V et al 2006 Phys. Rev. C 74 044602 

[3] Morita K, Morimoto K, Kaji D et al 2004 J. Phys. Soc. Jpn 73 2593 

[4] Adamian G G, Antonenko N V, Scheid W 1998 Nucl. Phys. A 633 409 

[5] Aritomo Y, Wada T, Ohta M 1999 Phys. Rev. C 59 796 

[6] Zagrebaev V I 2001 Phys. Rev. C 64 034606 

[7] Liu Z H and Bao J D 2004 Chin. Phys. Lett. 21 1491 
Liu Z H and Bao J D 2005 Chin. Phys. Lett. 22 3044 

[8] Li W F, Wang N, Li J F et al 2003 Eur. Phys. Lett. 64 750 
Li W F, Wang N, Jia F et al 2006 J. Phys. G 32 1143 

[9] Feng Z Q, Jin G M, Fu F et al 2005 Chin. Phys. Lett. 22 846 

[10] Feng Z Q, Jin G M, Fu F et al 2006 Nucl. Phys. A 771 50 

[11] Reiter P, Khoo T L, Lauritsen T et al 2000 Phys. Rev. Lett. 84 3542 

[12] Adamian G G, Antonenko N V, Scheid W 2003 Phys. Rev. C 68 034601 

[13] Grange P, Li J Q and WeidenmuUer H A 1983 Phys. Rev. C 27 2063 

[14] Li J Q, Tang X T, Wolschin G 1981 Phys. Lett. B 105 107 

[15] MoUer P, Nix J R, Myers W D et al 1995 At. Data Nucl. Data Tables 59 185 

[16] Oganessian Yu Ts, Utyonkov V K, Lobanov Yu V et al 2004 Phys. Rev. C 70 064609 




Figure 1: Angular momentum dependence of the partial cross section in the process of capture (a), 
fusion (b) and evaporation residue in 3n channel (c) for the system ^^Ca+'^^^ Cm at incident cm. 
energies 201.91 MeV and 211.91 MeV respectively, and evaporation residue excitation functions 
in 2n-4n channels compared with experimental data (d) 
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Figure 2: Evaporation residue excitation functions of the cold fusion ^^S'e+^°^5i and the reaction 
^^Ca +^^^ Bk in 2n-4n channels to synthesize superheavy element 117. 




Figure 3: Isotopic dependence of maximal production cross sections in In and 3n channels (up- 
per panel) and corresponding optimal excitation energies (lower panel) to synthesize superheavy 
element 117 in cold fusion reactions based '^^^Bi and in ^^Ca induced reactions. 
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Figure 4: Comparisons of calculated evaporation residue excitation functions in 2n-4n channels 

for the reaction systems ^^S'c +^^^ Cm(solid line) and "^^Sc +^^^ Cm(dashed line), ^^V ^-^^"^ Pu, 
^^Mn +^^^ U and ^^Co +^^^ Th to produce superheavy element 117. 



Table 1: Comparisons of calculated optimal excitation energies, maximal evaporation residue cross 
sections and transmission, fusion, survival probabilities at J = in ln-4n evaporation channels. 



Reactions 


Channels 








W,uriJ = 0) 

' ' SUT \ " " / 


Cpp/pb 


i^Ca +247 Bk 


2n 


26 


0.62 


7.19x10"'^ 


2.09x10-^ 


0.61 


^&Ca +247 Bk 


3n 


37 


0.88 


5.95x10"^ 


1.51x10-5 


0.86 


*^Ca +2-17 Bk 


4n 


47 


0.97 


9.69x10"'^ 


7.00x10-^ 


0.48 


i^Ca +248 Bk 


3n 


36 


0.86 


1.47x10"*^ 


5.13x10-5 


2.19 


i^Ca +249 Bk 


3n 


36 


0.85 


1.24x10"^ 


3.54x10-5 


2.14 


45 5^ +246 ^^ 


2n 


33 


0.45 


1.07x10"^ 


3.43x10-5 


0.45 


45 +246 Cm 


3n 


46 


0.81 


7.06x10"^ 


4.77x10-'^ 


0.84 


i^Sc +246 Cm 


4n 


58 


0.96 


6.91x10^4 


4.31x10-^ 


0.96 


^^Sc +248 Cm 


2n 


33 


0.39 


2.24x10"^ 


4.07x10-5 


0.77 


45 +248 Cm 


3n 


45 


0.74 


1.23xl0~4 


8.95x10-'^ 


2.06 


45 +248 Cm 


4n 


57 


0.94 


1.27x10"^ 


7.99x10-^ 


2.88 


51 Y +244 


2n 


27 


0.42 


2.41x10"^ 


1.20x10-^ 


0.29 


51 Y +244 


3n 


41 


0.79 


4.13x10"^ 


4.13x10"^ 


0.60 


51^+244 p^ 


4n 


54 


0.96 


7.44x10"^ 


2.17x10-^ 


0.43 


55 Mn +238 U 


2n 


27 


0.48 


5.72 xlO"'' 


1.03x10-^ 


0.07 


^^Mn +238 U 


3n 


41 


0.76 


6.39x10"^ 


3.14x10-^ 


0.068 


^^Mn +238 u 


4n 


56 


0.94 


1.54x10"^ 


9.58x10"^ 


0.038 


59 Co +232 pf^ 


2n 


27 


0.51 


1.47x10"*^ 


6.77x10-4 


0.101 


59Co +232 p^ 


3n 


42 


0.80 


1.13x10"*^ 


1.49x10-^ 


0.051 


59(70 +232 Th 


4n 


58 


0.95 


2.07x10-5 


4.31x10-^ 


0.019 


76 Se +209 Pi 


In 


13.5 


0.85 


1.54x10-^ 


1.64x10-^ 


0.025 


77 Se +209 Bi 


In 


12 


0.84 


1.55x10^^ 


1.03x10-2 


0.18 


7SSe +209 Pi 


In 


13 


0.85 


9.26x10-10 


4.71x10-3 


0.044 


79 Se +209 Bi 


In 


11 


0.83 


5.31x10-10 


3.64x10-2 


0.24 


80Se +209 Pi 


In 


12.5 


0.89 


2.08x10-11 


1.29x10-2 


0.0048 


S2Se +209 Pi 


In 


12 


0.91 


4.67x10-11 


3.14x10-2 


0.015 



in 



